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Abstract 



We compute the soft and virtual NNLO QCD corrections to Higgs production through 
gluon-gluon fusion at hadron colliders. We present numerical results obtained at the 
LHC and at the Tevatron Run II. 
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Although the Higgs boson is a crucial ingredient of the Standard Model (SM), it has so far 
eluded experimental discovery. LEP results |I| imply a lower limit of > 113.5 GeV (at 95% 
CL) on the mass of the SM Higgs boson, favour a light Higgs boson {Mh ^ 200 GeV) and 
suggest the observation of signal events at Mh — 115 GeV. 

After the end of the LEP program, the Higgs search will be carried out at hadron colliders, 
where the dominant production mechanism is gluon-gluon fusion through a heavy-quark loop. At 
the LHC [0 , gg fusion overwhelms the other production channels in the case of a light Higgs boson, 
and it still provides ~ 50% of the total production rate at Mh — 1 TeV. At the Tevatron 0, gg 
fusion gives ~ 65% of the total cross section in the mass range Mh = 100-200 GeV, although this 
production mechanism leads to an important discovery mode only at Mh ^ 135 GeV, when the 
decay channel H — > WW* opens up (if Mh ^ 135 GeV, the decay channel H —> bb suffers from a 
huge QCD background, and the decay rate if —> 77 is too low to be observed). 
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Figure 1: Parton luminosities as a function of r = s/s at the LHC. The inset plot shows the 
luminosities at the Tevatron as a function of the partonic centre-of-mass energy 

The NLO QCD corrections to Higgs boson production through gg fusion give a large effect [^]. 
Since approximate evaluations of higher-order terms suggest that their effect can still be size- 
able, the computation of the NNLO corrections is certainly important. 

A first step in this direction has been performed by two groups [|], ^ , who have computed the 
soft and virtual contributions to the NNLO partonic cross section (T{gg H+X) in the large- mtop 
approximation. Our calculation was done by combining the recent result for the two-loop 
amplitude gg —* H with the soft factorization formulae for tree-level ||^ {gg Hgg, Hqq) and 
one-loop [1^ {gg — > Hg) amplitudes. The independent calculation of Ref. uses a different 
method, and the results for the partonic cross section fully agree. 

In Ref. , we have combined the NNLO partonic cross section with the recent MRST2000 set 



of parton distributions |Tl| , which includes (approximate) NNLO densities, thus providing a first 
consistent estimate of the NNLO QCD corrections to the hadronic cross section. In the following, 
we briefly summarize our main results. 
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Figure 2: K-factors at the LHC (above) and Tevatron Run II (below): exact NLO result, NLO-SV 
and NLO-SVC approximations (left), and NNLO-SV and NNLO-SVC approximations (right). 

The soft and virtual terms give the most important contributions to the partonic cross section 
when it is computed at a value of the partonic centre-of-mass energy y/I that is close to the 
Higgs boson mass Mh- The hadronic centre-of-mass energy ^/s is related to the partonic one as 
s = Ts = X1X2S, where Xi and X2 are the momentum fractions of the two partons that initiate the 
partonic subprocess. The density of these partons is controlled by the parton distributions of the 
colliding hadrons. Since the parton distributions are strongly suppressed at large values of the 
momentum fractions Xi, X2, s is typically much smaller than s and, therefore, the soft and virtual 
terms can approximate the complete result also when the hadronic energy ^/s is not very close 
to Mh- At fixed Mh, the quantitative reliability of the approximation depends on the value of 
■s/s and on the actual value of the parton luminosity. The gg and [q + q)g luminosities are shown 
in Fig. |T]. At the LHC, where we use the reference value Mh = 115 GeV, the parton luminosity 
decreases by about two orders of magnitude when a/^ increases from Mh to "iMn (i.e. when r 
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increases from 10~^ to 10~^). At the Tevatron (see the inset plot), the parton luminosity decreases 
by almost two orders of magnitude when y/§ increases from \/§ = Mh = 150 GeV to 300 GeV. 
This observation suggests that the soft and virtual terms can give a good approximation at the 
LHC and a (slightly) better approximation at the Tevatron Run II. 

In Ref. 0, we introduced two approximations, named soft-virtual (SV) and soft-virtual- 
coUinear (SVC), of the partonic cross section. In the SV approximation only the contributions of 
soft and virtual origin are taken into account. The SVC approximation extends the SV approx- 
imation by including the next-to-dominant contribution (which has a coUinear origin |^) in the 
expansion around the region where = Mh. 

In the following we present numerical results for the hadronic cross section evaluated up to 
LO, NLO and NNLO. At the NLO we compare the exact result with those of the SV and SVC 
approximations. At the NNLO, we add the exact NLO result to our SV and SVC approximations 
of the NNLO corrections. All the results are presented in term of K-factors, which are defined 
as the ratio of the hadronic cross section over its LO value. The latter is computed by fixing 
the factorization and renormalization scales fip and fiu at the central value fJ^F = fJ'R = Mh- As 
a reference, we use the parton distributions of the MRST2000 set, with parton densities and as 
evaluated at each corresponding order. 

At the LHC, we find that the results do not strongly depend on the choice of parton distri- 



butions: very similar results are obtained by using the CTEQ5 set [T^, whereas larger (~ 10%) 
deviations appear by using the GRV98 set |TB[. At the Tevatron Run II, the dependence on 
the parton distributions is stronger than at the LHC. The CTEQ5 set gives a LO cross section at 
Mh = 100 GeV {Mh = 200 GeV) that is ~ 10% (~ 30%) smaller than the value obtained by using 
the MRST2000 set. Of course, this dependence affects the K-factors through their normalization 
with respect to the LO cross section. 

The K-factors at the LHC and at the Tevatron are shown in Fig. |^. The bands are obtained 
by independently varying fip and fiR in the range Mh/^ < fiF,fJ'R < 2Mh- The results show a 
reduction of the scale dependence at NNLO (from about ±20% at full NLO to about ±10% and 
±15% at NNLO-SV and NNLO-SVC, respectively). 

From the plots on the left-hand side of Fig. 0, we see that at NLO the SV approximation 
tends to underestimate the exact result, whereas the SVC approximation slightly overestimates it. 
Therefore, it is reasonable to expect the exact NNLO K-factor to lie inside the corresponding SV 
and SVC bands. In the case of the production of a light Higgs boson {Mh = 100-200 GeV) at the 
LHC (Tevatron), this expectation corresponds to i^' ~ 2.2-2.4 {K ~ 3), i.e. to an enhancement 
of the NLO cross section by about 10-25% (50%). 

Note that, at fixed Mh, Higgs boson production at the Tevatron is closer to threshold than at 
the LHC. This implies (see Fig. ^) that the approximation of the fixed-order results in terms of 
the corresponding SV and SVC contributions works better at the Tevatron than at the LHC. This 
also implies that, order by order in perturbation theory, QCD radiative corrections are larger at 
the Tevatron than at the LHC. In particular, QCD contributions beyond NNLO, due to multiple 
soft-gluon emission, can still be numerically relevant to Higgs boson production at the Tevatron. 
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